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SUMMARY 

This Interim report consists of a manuscript, "Receiver Design for 
Satellite to Satellite Laser Ranging Instrument," and copies of two 
papers we co-authored, "Demonstration of High Sensitivity Laser 
Ranging System" and "Semiconductor Laser-Based Ranging Instrument 
for Earth Gravity Measurements." These two papers were presented at 
the conference Semiconductor Lasers, Advanced Devices and 
Applications, August 21-23, 1 995, Keystone Colorado. The manuscript 
is a draft in the preparation for publication, which summarizes the 
theory we developed on space-borne laser ranging instrument for 
gravity measurements. 
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1 Introduction 

Precision laser ranging between a host satellite and a co-orbital passive subsatellite 
can be used to measure the gravity field distribution of the earth. The relative 
velocity between the host satellite and the subsatellite is a function of the gravity field 
distributions below the satellites. Precision range measurement can be achieved with 
a laser ranging instrument on the host satellite and retroreflectors on the subsatellite. 
Because the divergence angle of the laser beam is much smaller than that of radars, 
the signal propagation losses in space for a laser ranging instrument is much lower 
than radars. NASA has studied the use of a laser ranging instrument in the Gravity 
and Magnetic Earth Surveyor (GAMES) mission [l] [2]. The technique can also be 
used to map the gravity distribution of the moon and other planets. This paper 
presents a receiver design and performance analysis of such a laser ranging systems. 
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The basic laser ranging instrument consists of a laser transmitter, a photodetector, 
and some signal processing circuitry. Here we consider the use of a semiconductor laser 
with sinusoidal intensity modulation. The reflected laser beam by the subsatellite is 
collected by the receiver telescope and then detected by a photodetector. The received 
signal is then compared with the transmitted signal to determine the amount of 
phase shift which in turn gives the distance between the two satellites. The higher 
the modulation frequency, the larger the phase shift in the reflected signal and the 
finer the ranging resolution. On the other hand, the range measurement is periodic 
and the period is equal to that of the modulating signal multiplied by the speed of 
light and then divided by two. One can modulate the laser with several tones at 
successively higher frequencies to achieve both fine ranging resolution and wide range 
coverage. The receiver contains a narrow band filter for each tone and there should be 
no interference between channels. These type of laser ranging systems are relatively 
easy to built since they are direct detection system. 

The two satellites considered for NASA’s GAMES mission was 100 to 300 km 
apart and the relative velocity is expected to be about 1 m/s on average and changes 
very slowly. It is not practical to use coherent detection scheme in this application 
because the coherent time of lasers available today are shorter than or comparable 
to the round trip propagation time between the two satellites (0.67ms/100km). The 
laser transmitter being considered here are AlGaAs laser diodes at about 800nm 
wavelength. The advantage of this type of laser diodes are the small size, high effi- 
ciency, and relatively high output power. Sinusoidal intensity modulation is preferred 
because the receiver signal processing are similar to those in conventional Doppler 
radars and the technologies are mature and readily available. 

The photodetectors which render the highest receiver sensitivity are photomulti- 
plier tubes (PMT) and silicon avalanche photodiodes (APD). Both devices can be 
operated in either photon counting mode or analog mode. In photon counting mode, 


2 



a photocurrent pulse generated by a single photon absorption is detected by a dis- 
criminator as a discrete event. The use of the discriminator helps to reduce the effect 
of the photodetector gain randomness as well as the circuit thermal noise. On the 
other hand, the receiver is narrow band and should reject most of the thermal and 
photodetector gain noises. Therefore we may just use these photodetectors in analog 
mode and avoid the time walk introduced by the discriminators itself. 


2 The Link Equation 

We first consider single tone sinusoidal laser intensity modulation. The transmitted 
laser signal power, P x (t), can be written as 


P x (t) = m x P 0 cos(ujQt) + P Q 


( 1 ) 


with 0 < mi < 1 the modulation index, P a the average transmitted laser power, and 
uj 0 the modulating frequency. The peak power of the laser is (1 + m x )P 0 . 

The laser beam propagates to the subsatellite and a small fraction is reflected back 
by the corner cube. The reflected laser beam is assumed to be diffraction limited with 
the full divergence angle 9 C — 2 x 1.22A/<£ C . with (f> c the diameter of the corner cube 
and A is the wavelength. Here we assume the subsatellite has one corner cube on the 
rear and points to the host satellite. This may be achieved by designing the shape of 
the subsatellite like a bullet so that its orientation can be kept to a certain degree by 
the drag of the upper atmosphere. In the absence of atmosphere such as around the 
moon, the subsatellite can be a sphere covered with a number of corner cubes. We 
will expand our analysis to include multiple corner cubes in a later section. 

The fraction of transmitted signal captured by the receiver as a function of the 
range can be written as 


A{z) 


<t>c \ 2 ( <j>r V _ Vr(p 4 c4> 2 r 

a x z) \9 c z) ^ T 5.95A 2 o? x z A 


( 2 ) 
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where z is the range between the two satellite, a x the transmitted laser beam diver- 
gence angle, 0 r the diameter of receiver telescope, and r/ r the transmission efficiency 
of the receiver optics. Notice that the received signal power is proportional to the 
fourth power of the size of the corner cube and inversely proportional to the fourth 
power of the distance. 

The received optical signal power can be written as 

Pr(t) = A[z(t)}P x (t - 2^) + P b (3) 

c 

with P b the average background illumination. Note the range is a function of time. For 
convenience, we choose the time origin equal to the start of the observation interval, 
i.e. 

z(t) = ZQ + vt (4) 

with zo the range to the subsatellite at the beginning of the observation interval and 
v the relative speed between the two satellites. It is this relative speed which we have 
to estimate in order to determine the gravity field. As mentioned earlier, v is small 
and changes very slowly over time. The observation time interval in (4) is assumed to 
be small that v can be approximated as piecewise constant. The propagation losses 
can also be approximated as unchanged over the observation interval. Substitute (4) 
and (1) into (3), 

P r {t) = A(z 0 )miP o cos[uj 0 t + 9 0 + u D t] + A(z 0 )P o + P b (5) 


where 


a 2z° 
c 

_ v 
u>d — 2u/’o — 
c 


(6) 

(7) 


The term up represents the Doppler shift due to the relative motion of the two 
satellites. 
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3 Receiver Design 

The anticipated relativity velocity between the two satellites is only about 1 m/s and 
changes very slowly. The Doppler shift can be written as two terms 

v Au 

u>d — 2cjo — h 2a^o — (8) 

c c 

with v and /\v the average and the fluctuation of the velocity. The first term in (8) 
is just a constant frequency shift. The second term carries the information about the 
gravity field and needs to be determined with high accuracy. As mentioned earlier, 
it is assumed that v ~ 1 m/s and 

, 

2 ljq — T <§C 27 r (9) 

c 

with T the observation time. The receiver primarily detects the phase of the received 
signal. The average Doppler shift may be obtained from successive phase measure- 
ments over a much longer integration time. 

Gagliardi and Karp [3] showed a maximum a posteriori (MAP) phase estimator for 
an ideal photodetector and a Gaussian phase distribution. It is given as the solution 
to the integral equation 

rT 1 Ov 1 » 

9 = a 2 x(t)tan[-uj 0 (l -I )t 4- -9\dt (10) 

Jo 2 c 2 

where x(t) is the photon counting process output from the photodetector and a 2 is 
the variance of the phase. The equation for the MAP phase estimator given in (10) 
may be depicted as a phase tracking loop called a tan-lock loop [3]. However, it is 
very difficult to implement a tan-lock loop at the frequencies we are interested in 
(> 1GHz). 

The receiver we are considering consists of a narrow bandpass filter and a fre- 
quency and phase detector. It is simple in both concept and implementation though 
not mathematically optimal. The bandwidth of the filter should be as small as pos- 
sible but sufficient to accommodate the Doppler shift. Asin radio frequency (RF) 
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receivers, we further down convert the frequency of the received signal to some in- 
termediate frequency (IF) before the phase detection. The narrow bandpass filter 
and the subsequent phase detector become much easier to build at lower frequencies. 
It also enables us to use digital signal processing technique for the filtering and the 
phase and frequency detection. Figure 1 shows the block diagram of the receiver. 
The IF frequency should be low enough for easy signal processing but high enough 
to minimize the effect of the 1/f noise of the subsequent circuit. This type of laser 
ranging receiver has been used by Payne [4] for precision ranging of stationary targets. 

There are many kinds of digital phase and frequency detectors. We consider here 
a relative simple phase estimator given by 


9 — tan 1 


Im[Y(uj p )\ ) 
Re[Y{u p )\ ] 


tan 1 


£n=o y n sin(uj p nAt) 

£n=oynC0s(co p nA£) 


( 11 ) 


where i/ n s are the digitized samples of the down converted signal, Y (u>) is the Fourier 
transform of the sampled signal, A t is the sampling interval, N is the number of 
samples in the observation interval, and u> p is the predicted frequency of the received 
signal given by 

aip — ujjf + ujq ( 12 ) 


with Up the predicted average Doppler shift and u>i F the IF frequency. 

This phase estimator is similar to the one used by Payne [4] except that the 
target in our case is moving constantly and the exact frequency of the received signal 
is unknown. Because the Doppler shift changes very slowly over time, the prediction 
of the frequency can be very close once the receiver is tracked on. Next, we derive 
the effect of a small unknown frequency on the phase estimator. 

The average received signal can be written as 


y Q (t) - I s cos[(uj if + c j D )t + $q 


(13) 


with a ) D the unknown Doppler shift and 9q a constant phase. The sampled signal 
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after the A/D converter over a measurement interval can be expressed as 


JV 


y{t) = ~ n At) = E 6(t-nAt) ] ■ y 0 (t) ■ w(t) 


( 14 ) 


71=0 


where y n — yo (nAt) and xu(t) is the windowing function given by 

' 1, ifO<t<T 


w(t) = 


0, otherwise 


(15) 


The Fourier transform of y(t) is 

V» = 


At 

where V denotes convolution. 
Since 


~ E s (u-A 

A# V At 


k=-o c 


* Y„(ui) , W( a) 


(16) 


Yq(uj) — — - ^e? e °6[uj — (wjp + wd)\ + e ^ e °6[u + {ujjf + c^d)]} (.17) 


and 


W(u) = T • ^f^ ~ iuT/2 = T ■ sinc{uT/2 ) ■ e“ JwT/2 


(18) 


the Fourier transform of the signal can be written as 

[ui — (uijf + up) — kx;)T 


Y (w) — — — — { e ]6 ° V svnc 
to I *=-oc 


+ e 


-J0O 


£ 

/:=— oo 


S271C 


(a; 4- (w//r + U>d) — kjr:)T 


g -ji(w-(u;/F+w D )-fcgL)f 


The sample rate is assumed to be fast enough that there is no aliasing between each 
terms of the summations in (19) and all the k 0 terms may be filtered out. The 
observation time is set to some multiple of the period of the estimated frequency, i.e. , 


UpT — 2At7t, k = integer 

Substituting u = u p into (19) and dropping all the k ^ 0 terms 

'2lu v T — A uT\ -i &ut 


(20) 


/ ^ 7iTI s 


( A wT\ i( &*r 
smc I - — - 1 Jy 2 


a \ ( \ 2WT o 1 

e ^ — *o) + S i n c f — ) e jl ~ — * o1 


) 


( 21 ) 
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where 


A u — uj p — ( ujp ■+■ uq) 


( 22 ) 


The phase estimator becomes 
9 = tan 


( sin( 0o ) 

| sine 

/A uT 

K 2 

) — sine 

f 2uj p T— A u;T ^ 

V 2 ) 

11 

\ cos( A ; T 9 0 )[ 

sine ^ 

A uT\ 
2 ) 

+ sine ( 

2w p T — AwT \ 

2 Jl 

1 


Since u p T = 2kn, sinc(u p T — AuT/2) = — sin( AuT/2) /(u p T - AuT/2) 

1 


8 = tan 


-i 


, //) AuT . Aw T 2w p — Aw 

tan(6o f i — 

2w p “ A U/ 


-1 


AuT . ( 

tan(9 0 — ) 1 - 


Au; 


u p - A u t 


— tan 

If the frequency of the received signal is known, i.e. , Ao; = 0, 

9 = tan~ 1 [tan(9 0 )} = 9 0 


(23) 


(24) 


(25) 


When the frequency of the received signal is unknown, the estimated phase contains 
an additional term which is equal to the frequency difference between the predicted 
and the actual times one half the observation time. The unknown frequency also 
introduces a bias term in the estimated phase. The relative error in the estimated 
phase due to this bias is approximately 

Au; Au; 

Zb = T— = 

u p — Au ujf 

Therefore, the IF frequency should be much greater than the amount of the unknown 
Doppler shift. It also seen from (34) that the observation interval has to satisfy 
AuT/2 <C 7r in order to avoid ambiguity in the phase measurement. 



4 Receiver Signal to Noise Ratio 

The signal output from an ideal photodetector can be modeled as a Poisson random 
point process with the average counting rate with the average photocurrent given by 
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(51 


WO = <&- + h (27) 

hf 

where q is the electron charge, G is the average photodetector gain, 77 is the quantum 
efficiency of the photodetector, hf is the photon energy, and I d is the photodetector 
dark current. The amplitude of the sinusoidal signal is obtained by substituting (5) 
into (27), as 

I. = qG^-A(z 0 )P o m 1 (28) 

nj 

The total noise output from the photodetector consists of the shot noise associ- 
ated with the photon detection and the circuit noise. The latter can be modeled as 
independent and additive Gaussian noise. The former is signal dependent and the 
exact distribution function is not readily available. However, as shown in [6], the pho- 
todetector shot noise can be considered as independent Gaussian noise if the receiver 
is narrow banded and the signal to noise ratio is high. The one sided noise current 
spectral density can be written as [7] 

N n = 2qG 2 F[<I pd (t)>+I d \ + N amp 

= 2qG 2 F[^A(z 0 )P o + ^P b + I d ] + N amp (29) 

where F is the excess noise facto, N amp is the noise density of the preamplifier, and 
< Ipd{t) > is the average DC photocurrent. 

If a photomultiplier tube (PMT) is used, N amp can often be neglected because the 
average photodetector gain is large. The excess noise factor of a PMT is given by 
Fpmt ~ m/(m — 1) with M the average gain of each dynode (m ~ 3 — 5). If an APD 
is used, the excess noise factor is given by [8] Fapd — KfjG -I- (2 — l/G)(l — kef f) 
with k e fj the ratio of the hole to electron ionization coefficients. The circuit noise has 
to be considered in this case. The circuit noise consists of primarily the preamplifier 
noise, which is often given by the manufacture or can be measured directly. 
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The signal to noise ratio (SNR) at the output of the narrow band filter is given as 


SNR = 


N n B n 

\ [hf A ( Z o) p o^f B 

2F[%A(z 0 )P o +j f Pi, + *}}B n + %& ' 


with B n the one sided noise bandwidth of the filter. 


(30) 
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